Zener tunneling in the electrical transport of quasi-metallic carbon nanotubes 
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We study theoretically the impact of Zener tunneling on the charge-transport properties of quasi- 
metallic (Qm) carbon nanotubes (characterized by forbidden band gaps of few tens of meV). We also 
analyze the interplay between Zener tunneling and elastic scattering on defects. To this purpose we use a 
model based on the master equation for the density matrix, that takes into account the inter-band Zener 
transitions induced by the electric field (a quantum mechanical effect), the electron-defect scattering 
and the electron-phonon scattering. In presence of Zener tunnelling the Qm tubes support an electrical 
current even when the Fermi energy lies in the forbidden band gap. In absence of elastic scattering (in 
high quality samples), the small size of the band gap of Qm tubes enables Zener tunnelling for realistic 
values of the the electric field (above ~ 1 V//^m). The presence of a strong elastic scattering (in low 
quality samples) further decreases the values of the field required to observe Zener tunnelling. Indeed, for 
elastic-scattering lengths of the order of 50 nm, Zener tunnelling affects the current/voltage characteristic 
already in the linear regime. In other words, in quasi-metallic tubes, Zener tunneling is made more visible 
by defects. 



I. INTRODUCTION 

Single-wall carbon nanotubes (SWNTs) are quasi one- 
dimensional w/ires of great interest for future electronic- 
devices applications. A SWNT can be constructed by rolling 
up a graphene sheet and its geometry is univocally speci- 
fied by a pair of chiral indexes (n, m)-'-. The relationship 
between n and m defines three groups of SWNTS: arm- 
chair {n,m = n), zigzag {n = or m = 0) and chiral 
nanotubes. The indexes {n,m) also determine the SWNT 
electronic structure. Ifm — rtisnota multiple of 3, SWNTs 
of diameters in the 1-2 nm range are semiconductors with 
a band gap A larger than 0.5 eV. Only armchair SWNTs 
were believed to be truly metallic. Indeed, because of the 
curvature, zigzag and chiral SWNTs with m — n multiple 
of 3, present a small band gap of few tens of meV— , and 
for this reason are called quasi-metallic (Qm)^. However, 
recent transport measurements on chirality-identified tubes 
showed that also armchair tubes are Qm with A ~ 10 to 
100 meV, likely because of electrons correlation^.. 

The description and interpretation of the transport prop- 
erties of SWNT devices, are often based on a semi- 
classical approach such as the Boltzmann transport equa- 
tion (BTE)^. However, the BTE can not describe some 
cases such us strong localization regimes^. Another effect 
that can not be described by the Boltzmann transport ap- 
proach is the Zener tunneling (ZT)^, which is the electric- 
field-induced tunneling of carriers from one band to another 
through the forbidden energy gap A. In absence of ZT, if 
the temperature is smaller than the gap (fc^T ^ A), only 
one type of carrier is responsible for charge transport: when 
the Fermi level, ep, is in the conduction (valence) band, the 
current is carried by electrons (holes). If ep is in the gap 
region charge-transport does not occur. On the contrary, 
in presence of ZT, the current is carried simultaneously by 
both electrons and holes and a non-zero current exists even 
if ep is in the gap region. 

The critical electric field for which the ZT affects 
transport properties, can be estimated considering the inter- 
band tunneling probability across the gap of the SWNT 



hyperbolic bands. In absence of scattering, the transition 
probability as a function of the source-drain electric field 
EsD is given by = ex-p{—TrEz/EsD}^ where 
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and vp ^ 10^ cm.s"^ is the Fermi velocity of graphene 
(and of Qm nanotubes far from the gap). In semiconduct- 
ing tubes (A ~ 500 meV), E^ ~ 95 V//im is unrealistically 
large. At experimentally-accessible electric fields, ZT is ab- 
sent, and in a diffusive regime a semi-classical picture is 
sufficient to describe transport properties. Instead, in Qm 
tubes (A ~ 50 meV) E^ ^ 0.95 \//fim, a value that nan- 
otubes can sustain^. In this case ZT could be observed. 

In this work, we study theoretically the current voltage 
characteristics, of Qm nanotubes both in the linear and high 
field regime and we analyze the interplay between Zener 
tunneling and defect (or impurity) scattering. We consider 
long Qm tubes with an energy band gap A = 60 mcV (sim- 
ilar to the values measured in experiments^;!) in a diffusive 
regime, i.e. when the tube length is much longer than the 
defect- induced elastic scattering lengths {1^) of 50 nm and 
300 nm presently used. The approach is based on the solu- 
tion of the master equation for the evolution of the density 
matrix. This approach correctly describes inter-band Zener 
tunneling. 



II. GENERAL ASSUMPTIONS 

We take into account the elastic scattering of electrons 
with defects (we consider scattering of electrons by short- 
range potential) and the inelastic scattering with optical 
phonons. We do not consider the possible occurence of 
weak or strong localization effects^. Elastic and phonon 
scattering are assumed to be incoherent and, thus, the 
phase coherence length < l^. Scattering with acous- 
tic phonons plays a minor role since it is characterized by 
scattering lengths that are much larger than the le used 
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here (see e.g. Ref. [^. We suppose that the tubes are sup- 
ported and that, for the regime under study, the phonon 
population is equilibrated at the temperature of the sub- 
strate. We assume that the doping and the potential drop 
is uniform along the channel, and that the electron popula- 
tion of the valence and conduction bands is independent on 
the position in the channel. Namely, we consider the bulk 
transport properties of an infinite SWNT under a spatially 
homogeneous electric-field applied along the tube axis. 

We remind that at high bias conditions in small gap nan- 
otubes, the scattering with optical phonons can lead to a 
current saturation^ and to nonequilibrium optical-phonon 
population (hot phonons}i2iii. In suspended metallic tubes, 
this can lead to the observation of negative differential re- 
sistanc o^i-*-^ . This kind of phenomena are not included in 
the present model and their inclusion would lead to increase 
of the differential resistance for high electric field (this point 
will be discussed in Sects. \\\/[ IV Bp . 

Electronic transport measurements in Qm nanotubes 
were reported in The two papers deal with trans- 

port regimes which are quite differents from the present 
one. Because of this, the results of 01^ can not be directly 
compared with the present ones. In particular, in both 
the elastic scattering lengths are much longer than in the 
present work and, in both cases, they are comparable with 
the length of the tubes. Because of this, Ref. 4 can observe 
interference effects between the contacts and Ref. can 
interpret the results by assuming an electronic distribution 
non uniform along the length of the channel. On the con- 
trary, in the present work we are in a diffusive regime: the 
scattering length If, ^ L, where L is the length of the tube, 
and Lcf, ^ L, where is the phase coherence length. In 
this diffusive regime, we can assume a uniform distribution 
of the electrons along the tube and that interference effects 
between the contacts are not observable. 



III. DENSITY MATRIX IN AN UNIFORM ELECTRIC 
FIELD 

A. Hamiltonian 

The graphene Hamiltonian for the 7T-bands close to the 



K point (Dirac point) of the Brillouin zone (BZ) 

^■20 = vpo- ■ p 



,-13. 



(2) 



where vp is the Fermi velocity, a = {<Tx,(7y) with and a-y 
the Pauli matrices, p is the momentum operator in the two- 
dimensional BZ of graphene. 7^2D is a 2x2 matrix which 
describes the tt and tt* bands. Starting from Eq.(l2]), one 
can derive an effective Hamiltonian H for a SWNT: 



(3) 



Px is the momentum operator in the one-dimensional BZ of 
the nanotube (a momentum parallel to the tube axis) and 
ko is a constant fixed by the energy Gap A 2hvFko. The 



eigenstates of H are Bloch states and have the form 



(4) 



Here L is the length of the system and \Uk.a). is the periodic 
part of the Bloch state: 



|f/fe,a> 
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where 6 = arctan(fco/fc)- a = 1 indicates positive energy 
(conduction band) whereas a ~ —1 describes states with 
negative energy (valence band) with respect to the Dirac 
point. The eigenvalues of H read: 



(6) 



We now apply a spatially homogeneous and static (source- 
drain) electric field Esd along the tube axis. The Hamilto- 
nian, in the vector-potential gauge, becomes : 



y. = vf[(Jx{Px + eA{t)/c) + cjyhko] 



(7) 



Where A{t) = —cEsot is the vector potential component 
along the tube axis, and e = |e| the electron charge. The 
instantaneous eigenstates of the time dependent Hamilto- 
nian, Eq.d?]), are given byii : 



(8) 



\k,a} =e-''^^''>^^^^\k,a) 
with eigenvalues ea{k). 

B. High field: non-linear regime 



A framework in which quantum effects such as Zener tun- 
neling appear naturally is the Liouville-Von Neumann Mas- 
ter Equation (ME)i^. Within this formalism our variable is 
the one-electron density-matrix operator p. The time evo- 
lution of p in absence of scattering processes is given by: 



dp 

dt 



■m,p]. 



(9) 



Since we are dealing with a spatially homogeneous elec- 
tric field applied to an infinite nanotube, only the elements 
of the density matrix diagonal with respect to k will be 
relevant—. We use the instantaneous eigenstates (Eq.([5|)) 
as basis. Then, from Eq.jQ]) one finds the time evolution of 
Pa!3{k,t) = {k,a\p\k,l3}^: 



dpai3{k,t) eEsD dpap{k,t) Aeal3 



dt 
eEsD 
ih 



dk 



ih 



Pal3{k,t) 



'^{Ran{k)pni3{k,t) - Par,{k,t)R,^fi{k)) 



(10) 



where Ae^^ = ea{k) — e/3(fc) and i?^,, is defined by: 

i?a^(fc) = (f/fc,a|i^|f/M>(l - '^ar,). (H) 
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The Ran term introduces the inter-band transitions. Using 
perturbation theory for the derivative with respect to k we 
can show that: 



{k,a\v\k,r]) 

where v is the velocity operator defined by: 

1 

For our system takes the form: 

1 k 



2 fc2 + fc2 



(1 - Sajj). 



(12) 



(13) 



(14) 



The full evolution equation is obtained by adding the 
collision integral I[p]afs to the right-hand side of Ea. (|10p . 
The collision integral results from the sum of two indepen- 
dent contributions due to inelastic electron-phonon scatter- 
ing and elastic electron-defect scattering. Electron-phonon 
scattering is treated as in Ref.-*^^ and read: 



e—ph 

a/3 ~ 



k' 



2L 



k' rjv 

y.paf}{k){l - Pr,r,{k')) 



(15) 



where W^^-^^'r^' electron-phonon scattering probability 
is evaluated using Fermi golden rulei^: 

W^fc',->fca = yl.9,";a(fc',fc)PK + 1/2 ± 1/2) ^^g^ 

x<5(e„(fc)-e^(fc')±K) 

with gj^^ the coupling matrix elements and v the phonon 
branch, fvjj'^^ is the phonon energy with q =\k' — k\ and 
is the phonon population. The upper (lower) sign in the 
delta function corresponds to phonon emission (absorption). 

The elastic scattering by defects is supposed to be inco- 
herent. In this case the collision integrals for the diagonal 
and non-diagonal elements [a ^ (3) of the density matrix 
are given by: 



k' ci^ka 



[Paa(fc')-Paa(fc)] (17) 



1[pU = E(^fc'"->fe- + W^}p^^p)p^p{k). (18) 



The transition probability, within the Fermi golden rule, is 
defined by : 

= '^-^^^\{k\a\VD\k,a)\H{ea{k') - £„(fc)) 



(19) 



Where No is the density of defects (that is the number of 
defects per unit length) and Vd is the scattering potential. 

Finally, the master equation, that describes the time evo- 
lution reads: 



dpapjk) ^ eEsD dpap{k) _ Ae^ff 
h dk ih 
eEsD 



dt 



Papik) 



+ 



ih 



Y,{Rcv{k)p.Aky (20) 



e—ph 



Pc^{k)Rn0{k))+l[p]%+l[pr-^ 

For a given value of the electric field Esd, a steady-state 
solution of Ea. (|20p is reached for t oo. In this work, the 
steady-solution is obtained by numerical integration in time 
using a finite-difference method. 

The electronic population (that is the diagonal elements 
Paa) determines the density of electrons N (number of elec- 
trons in excess per unit length of the tube) by: 



N 



9s9v 
L 



E 



(21) 



where g,, = 2 and = 2 are respectively the spin and valley 
degeneracy. is a constant of motion of Eg. ([20]) . We will 
show simulations done at different values of the Fermi level 
ep, that is by fixing the initial value of to: 



N 



9s9v 
L 



E 



E/(^"(^)"^f^)-i 



(22) 



Where /(e) is the Fermi-Dirac function, computed here for 
T = 50 K. 

The charge current is given by : 



1 = 



-9s9ve 



^^Vafi{k)pap{k), 

k a0 



(23) 



where Vap is defined by 



v^pik) = {k,a\-[H,x]\k,P) = (fc,a|-[H,x]|fc,/3) (24) 
n n 

Using the definition of the velocity operator given in 
Ea. ()13p . one can show that 



hVapik) 



del{k)ldk 



if a = /3 



-iRai3{k)Aeaf3 \f a ^ P 



(25) 



Finally the steady current will be evaluated for the steady- 
state solution. 

An interesting aspect of the master equation, is that 
we can recover the standard Boltzmann transport equation 
(BTE), i.e. the semiclassical description, setting i?,,^ of 
Ea. ()20p to zero. In this case the diagonal elements of the 
density matrix evolve according to: 



dpaa{k) eEsD dpaa{k) 



e—ph 



(26) 



dt h dk 

and the off-diagonal elements of the steady solution vanish 
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C. Linear regime 



D. Zero-field electrical conductivity 



In this section, we derive a set of simplified equations 
(which can be evaluated analytically) valid for small values 
of the electric field Esd within the linear regime. In the 
linear regime the applied field will cause a small deviation 
of the equilibrium distribution {Esd = 0). The steady state 
density matrix (with Esd 7^ 0) can be written as: 



Papik) 



(27) 



where pf^ is the density matrix element in absence of field 



{Esd ~ 0) and p''^^ the deviation linear in Esd- P^'g is 
the equilibrium distribution which is diagonal and equal to 
pi"^ = Sa0f{ea{k)~'€F)- Inserting Eq.^ in Eq.^. and 
keeping all the terms linear in Esd, one obtains the solution 
for the diagonal elements of the density matrix: 



(0) 



pi'lik)=Trik) 



eEsD dp^^lik) 



h dk 

with T,., the momentum relaxation time, defined byiS; 

^ TJ2<^^k'Mk,k') 



Tr{k) L ^ 



(28) 



(29) 



where ip{k,k') — 2 for back-scattering processes and for 
forward scattering ones. Indeed, after the scattering with a 
defect (or a phonon), the electron can either maintain its 
propagation direction, which is forward scattering (fs), or 
reverse it, which is back-scattering (bs). The off-diagonal 
elements (a ^ (3) are given by: 



zEsDRgpjk) 



Aeap - ih^tot(pal{k) - p^°j(fc)) 



where "ftot is defined by : 

7tot(fc) = lei{k) + 7,%(fc) + 7,^pjfc) (31) 

where 7pj, 7^^,^ and 7gp;j, are respectively the elastic, hole- 
phonon and electron-phonon scattering rate for the equilib- 
rium distribution and defined by: 



lel{ 



ka^k'a + ^kp^k' p) 



(32) 



lephik) = ^Y.T.(^l^'V^ka+Wi:,,_kp)pi°]{k') (33) 



k' rjV 



(34) 



k' r^u 



The zero-field conductivity, a" = Muiesd^o I / Esd< can 
be obtained by inserting the steady state solutions of the 
linear regime, Eas. (f28|) and (|30p in the current expression, 
Eq.dig): 

a" = a",+a: (35) 

where cr^ is the conductivity obtained from the diagonal 
elements of the density matrix. It is given by: 



k a 



Tr{k) ( dp^Sl{k) 



dk 



(36) 



(T^ is the semi-classical conductivity of the Boltzmann equa- 
tion and at zero temperature it reduces to: 

o e^vl^[kF)_ 



^(fcF)D(ei.) 



(37) 



with hkp the Fermi momentum and D(e) the one- 
dimensional density of states (DOS) per unit length. For 
the lowest sub-band of the Qm nanotube^S: 



- (A/2)^ 



= e(|e|-A/2) 



(38) 



where Dn 



is the constant DOS (per unit length) 
of the lowest sub-band of a metallic tube. <d{x) is the step 
function which equals 1 for a: > and otherwise. In the 
case of a metallic tube the conductivity, Ea.(|37p. can be 
written as: 



<J°n = -^leikp) 



(39) 



(30) where leikp), the elastic scattering length, is given by: 



leikp) = VpTrikp). 



(40) 



The current due to Zener inter-band is given by a°, 
that arises from the off-diagonal part of the density matrix, 
Eq. 

" = e„(fc)-e,(fc) 



L 



k aP 

f^ltotik) 



[eo.{k)-ep{k)f + [h^tot{k)] 



;{l~5^0). (41) 



Notice that this expression is similar to the off-diagonal 
Kubo formula for the optical conductivity^-"^ but with a finite 
value of the parameter ^tot- When ^tot = 0, the Zener con- 
tribution {(Tz) vanishes and the total conductivity is equal to 
that of the semiclassical Boltzmann approach (cr^). Thus, 
in the linear regime, the Zener tunnelling contributes to the 
transport only in presence of scattering events that limit 
the carrier lifetime. On the contrary, at finite field, Zener 
tunneling can occur in the ballistic case (see formula right 
above Eq.©). 
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FIG. 1: Electronic band structures of a quasi-metallic (A 7^ 0) 
and a metallic carbon nanotube (A = 0). For the Qm tube, 
we distinguish between valence and conduction bands [labeled 
(-1) and (-1-1)]. On the contrary, for the metallic tube it is 
more convenient to distinguish between bands in which the 
sign of the electron velocity is the same, that is negative (n) 
or positive (p). 



IV. TRANSPORT CALCULATIONS 

We performed calculations for a metallic tube and a quasi- 
metallic tube (Qm) with an energy band gap A = 60 mcV 
(similar to the values measured in experiments^'^). The 
tubes are subject to a spatially uniform electric field and the 
number of carriers (electrons or holes) is varied by changing 
the Fermi level ep. In the metallic tube the dispersion rela- 
tion is linear e±(fc) = zthiipk, with vp ~ 8.39 x 10' cm/s, 
in a Qm tube the bands are hyperbolas, Eq.([5]). 

The transport properties of the Qm nanotube are explored 
by a semi-classical and a quantum-mechanical approach. In 
the semi-classical case, we do not take into account the 
inter-band transitions due to the electric field by setting to 
zero the off-diagonal elements of the density matrix {pai3 
with a /3 ). The evolution of the system is then governed 
by the BTE, Eq. ([26|1. In the quantum approach, we take 
into account the inter-band transitions due to the electric 
field by considering the off-diagonal elements of the density 
matrix. The evolution is then governed by the ME, Eq. ()20p . 

For the metallic tube, the ME that governs the quantum 
evolution can be obtained in the limit A ^ 0. In this 
case, the coupling [(Eq. (fTTjl]. that promotes Zener 
transitions, diverges at /c = 0, and the electrons that cross 
the fc = point tunnel with 100% probability from the 
valence (conduction) to the conduction (valence) band. As 
a consequence, for metallic tubes, the quantum-mechanical 
approach (ME) reduces to the semiclassical approach (BTE) 
if we change the meaning of the band indexes in Eq. [26l 
In particular, while for Qm tubes we distinguish between 
"valence" and "conduction" bands, for metallic tubes we 
distinguish between bands with positive or negative velocity 
{de/dk < or de/dk > 0), see FigH 

We consider short-range scatterers and the electron- 



defect coupling alVolfc, a) I is taken to be a constant 
and /c-independent. Under this hypothesis in a metallic 
tube the elastic scattering mean free path = t"'vp is 
^-independent. The value of alyol^', a)| is determined 
fixing I'J" = 50 nm and ;™ = 300 nm. In the Qm tubes, 
|(fc', ajVolfc, a)| is also assumed to be a constant equal to 
that of the metallic tubes. Note, however, that, in the Qm 
case, the scattering mean free path become fc-dependent. 
Indeed, because the presence of a gap, the density of states, 
that determines the scattering rate, is energy dependent. 

Only three optical phonons can play an important role in 
transport phenomena at high bias: two (longitudinal and 
transverse) correspond to the graphene optical phonons at 
r with symmetry E2g, and one corresponds to the graphene 
phonon at K with symmetry A'j^^ -'^^-'^-'^'^^. The energies of 
these phonons are respectively huj^ = 200 mcV and huj^ ~ 
150 meV. For our calculation the electron-phonon coupling 
matrix elements g'j'faik' , k) ~ 5^ are taken constant and in- 
dependent of k and k' , with A — bs, fs. The label 6s stands 
for a back-scattering process (f,,,, (A;)waa(fc') < 0), and fs 
stands for a forward scattering process {v,p^{k)vaa{k') > 0). 
In a metallic tube, the corresponding scattering lengths, 
Z^, are proportional to the tube diameter (P-^, that we as- 
sume to be 2 nm. In this case, = l^^ = 451.38 nm, 
= 183.74 nm and = 0. For the Qm tubes we use 
the obtained in the metallic case. The phonon popula- 
tion is set to zero (cold phonons, only phonon emission is 
treated). We do not consider the generation of hot phonons, 
that occurs at high bias in nanotubes, when the elastic scat- 
tering length is much larger than the phonon one, . In 
the present case, this is justified by the strong elastic scat- 
tering used in the present model that competes with the 
phonon scattering, as in graphene^^. 



V. RESULTS AND DISCUSSION 



We now show the results obtained for two different val- 
ues of the elastic scattering length, namely 1^ = 50 nm and 
le = 300 nm. The goal is to understand in which condi- 
tions the Zener tunneling is experimentally observable. We 
will compare the results obtained with the semi-classical (or 
Boltzmann) approach with those obtained with the quan- 
tum treatment (master equation). In the semi-classical ap- 
proach the Zener tunneling is not allowed. Thus, the semi- 
classical and quantum results will differ when Zener tunnel- 
ing is relevant. Results for metallic tubes (for which the 
distinction quantum versus, semi-classical does not hold) 
are shown as a comparison. The results shown for a finite 
value of the source-drain electric field, Esd, are obtained 
from numerical solution of the transport equations. The 
results for Esd ~ are derived analytically in Sec. INI Dl 
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FIG. 2: (Color online) Finite field conductivity a = I/Esd 
versus the Fermi energy ep for difTerent values of the field 
EsD- The two panels report results obtained with two differ- 
ent values of the elastic scattering length l^. aZi and a"° are 
the zero field conductivity of the metallic tube and of the Qm 
tube, respectively, cr^ is the component of cr° due to Zener 
tunneling. 



A. Conductivity and current as a function of the Fermi 
energy 

Fig. [2] reports the conductivity a of the Qm tube as a 
function of the Fermi level ep for various values of the elec- 
tric field EsD- In Fig. [2] is the zero-bias {Esd = 0) 
conductivity calculated in the linear regime for the temper- 
ature T = 50 K. (7° is the component of (t° due to Zener 
tunneling (cr" is a well defined quantity, see Sec. INI D|) . and 
(tJ^ is the zero-bias conductivity of a metallic tube. 

For the metallic tube, the conductivity (t°„ has a finite 
value which is proportional to the elastic scattering length 
If. and which does not depend on ep (this is true as soon 
as ep does not meet the higher subbands of the SWNT, 
which are not considered in this study). For the Qm tube, 
the conductivity tends to be equal to the metallic one for 
cp ^ A/2: when the doping is sufficiently high there is only 
one type of carriers, the details of the electronic bands near 
the Dirac point are not relevant, the transport is properly 



described by the semi-classical Boltzmann approach and the 
Qm and metallic tube behavior are indistinguishable. On the 
other hand, for ep < A/2 the presence of the electronic gap 
A acts as a barrier and the Qm conductivity suddenly di- 
minishes. In particular, from Fig. [2] for ep < A/2, a'^ ~ cr^ 
meaning that the Qm conductivity is entirely due to Zener 
tunneling (inter-band) and that the Boltzmann conductiv- 
ity (intra-band) vanishes. Fig. [2] also reports the finite field 
conductivity a = I /Esd for different values of the electric 
field. Also at finite field, the conductivity is entirely due to 
Zener tunneling for small values of ep. 

The comparison of the results obtained for different values 
of the elastic scattering lengths, 1^, is very remarkable. In- 
deed, for small doping [ep < A/2), the zero-bias conductiv- 
ity for the smallest scattering length (cr'^ ~ 0.82 S.m^^ for 
le = 50 nm in Fig. [2]), is higher than the one for the longest 
scattering length, ((t° ~ 0.2 S.m^^ for 1^ 300 nm). More- 
over, most important, the ratio between the high doping (t° 
[ep >0.2 eV) and the small doping cr" is much higher for 
the smallest scattering length than for the longest one. 

This counterintuitive behavior can be understood by 
considering the electronic scattering rate (or broadening) 
ltot{k), defined in Eq. [3TJ jtot{k) represents the uncer- 
tainty with which the energy of an electronic band can be 
defined. Fig. [3] reports the electronic bands for the Qm 
tube: the bands are represented as lines with a finite thick- 
ness (gray area in Fig. [3]) equal to jtotik). The increased 
value of jtotik) in the vicinity of the gap is due to the elas- 
tic scattering. The broadened electronic bands of Fig. [3] 
overlap near k — for 1^ — 50 nm (low quality samples), 
meaning that, because of energy-uncertainty, the electrons 
do not see the gap and can tunnel from a band to the 
other even if the applied field is small. On the contrary, for 
le = 300 nm (better quality samples) the bands do not over- 
lap and for small electric field the Zener tunneling can not 
take place. In Fig. [31 the sudden increase of the broaden- 
ing for |fc| > 0.03 A^^ is due to the scattering with optical 
phonons. This scattering mechanism is active only when 
|e(fc)| > A/2 + huj^ and |e(fc)| > A/2 + hw^ , which is far 
from the gap. Zener tunneling is, thus, activated only by the 
elastic scattering (with defects) and not by the scattering 
with optical phonons. 

Finally, to see the things from a different perspective. 
Fig. |4] reports the current as a function of ep for various 
values of the electric field Esd- For large doping (large 
ep), the Qm tube current calculated with the semiclassical 
and with the quantum approaches tend to be equal. Also, 
both currents tend to be equal to the metallic tube one. 
The semi-classical approach predicts a zero current in the 
gap and underestimates the current in a range of ep that 
increases with Esd- In the region of the gap, the current 
is thus entirely due to Zener tunneling. 



B. Current as a function of the electric field 

To analyze the results as a function of the electric field 
Esd, we consider two characteristic fields: E^ and Eph. E^ 
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FIG. 3: (Color online) Broadening of the electronic bands 
of the Qm tube for different values of the elastic scattering 
lengths le- The thin (red) solid lines are the energy bands 
Efe. The gray area is the region delimited by the curves ± 
hjtot{k) /2, where "ttot(k) is the electronic scattering rate. For 
le ~ 50 nm the grey areas of the two electronic bands are 
superimposed in a region highlighted in black. 



FIG. 4: (Color online) Current versus ep for different values 
of the electric field with le = 300 nm and 1^ = 50 nm. Calcu- 
lations for the Qm tubes are done with the semi-classical or 
with the quantum approach. Calculations for metallic tubes 
are also shown. Note that, for ep — ^ 0, the current vanishes 
in the semi-classical approach, while it has a finite value in 
the quantum approach. 



is defined in Eq. ([T]) and it is the value of the electric field for 
which, in the absence of scattering, the Zener scattering is 
expected to be sizable. Eph is defined as Eph = huj^/{elg) 
and it is the value of the field for which optical-phonon emis- 
sion become sizable compared to elastic scattering. Eph is 
higher for smaller values of l^- To understand this defini- 
tion, one should consider that optical phonons have a finite 
energy hu'^ . In order to emit an optical phonon an elec- 
tronic carrier needs to attain a minimum activation energy 
corresponding to hoj'^ . This is possible if the carrier has ac- 
celerated in the electric field Esd at least through a length 
lact — I [eEsD^- The process is dominated by the 
zone boundary phonons since lo^ < uj^ . For small electric 
fields, Esd ^ Eph, 1^ ^ lact and the scattering with de- 
fects prevents the carriers to accelerate and, thus, to reach 
the activation energy. On the other hand, if Esd S> Eph, 
le 3> lact and the carriers can attain the activation energy 
before scattering with the defects. Thus, for small electric 
field, elastic scattering is the dominant scattering mech- 



anism (in this regime the current-voltage curve of metallic 
tubes is linear-). For electric fields higher that Eph the scat- 
tering with optical phonons is activated. This is associated 
with a decrease of the differential conductivity (sub-linear 
current-voltage curve) which can lead to a saturation of the 
current at high drain-source voltages in metallic tubes with 
long elastic scattering lengths^. The scattering with optical 
phonons can lead to an anomalous increase of the optical- 
phonon population (hot-phonons)i2iii which, in turn, will 
further diminish the high-bias conductivity. 

Fig. [5] shows the current as a function of the electric 
field for different values of ep. At high doping, the results 
for the semi-classical (BTE) and of the quantum (ME) ap- 
proaches are very similar and approach the metallic results. 
The two approaches provide significantly different results for 
sufficiently high electric fields and small values of the Fermi 
level. In particular, for = the current is exactly zero 
within the semi-classical model, while the current has a fi- 
nite value, which increases with the electric field, within the 
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FIG. 5: (Color online) Current versus, electric field Esd 
for le — 300 nm (top panels), le — 50 nm (bottom panels). 
Thin continuous lines (red): quasi metallic tube by varying 
ep from to 0.5 eV with 0.05 eV step. Thick dashed lines: 
metallic tube. Calculations are done with the semi-classical 
(left panels) or the quantum (right panels) approaches. 
and Eph are respectively the field from which Zener tunneling 
and the optical phonon scattering become important. E^ — 
1.37V//im, Eph ~ 0.53 V/^m for = 300 nm, Eph ~ 3.2 
Y/nm for le = 50 nm. Note that for ep — 0, the current in 
the semi-classical model is exactly zero. 



quantum approach, Fig. [5l For ep = and ep = 0.05 eV, 
within the quantum approach the current versus, electric 
field is a super-linear curve for Esd ~ E^- Such a be- 
havior was observed in graphene^* and in large diameter 
nanotubes-^=, and identified as due to Zener tunneling. 

Let us compare the results for the different scattering 
lengths le (Fig. [5]). Zener tunneling should be more visible 
for le = 50 nm, than for Ip = 300 nm, for two reasons. 
First, for 1^ = 50 nm, the Zener current is already rele- 
vant for electric fields smaller than E^ (this is explained in 
the discussion related to Fig. [3]). Second, for 1^ = 50 nm 
Eph > Ez, while for le = 300 nm Eph < E,. Now, the 
scattering with optical phonons is associated to a sub-linear 
current which tend to mask the presence of the super-linear 
Zener current. For le = 50 nm, Eph > E. means that this 
"masking" does not occur in the region of interest. For le ~ 



300 nm, Eph < means that this "masking" is already 
active in the region where Zener is active. Beside this, one 
should also consider that in the present simulations we are 
not including the possibility of hot phonons. The inclusion 
of this effect would decrease the differential conductance 
for Esd > Eph- This should not change the le = 50 nm 
curves but should further mask the Zener current in the le = 
300 nm case. 

Finally, in by Fig[6]and[7]we present the two-dimensional 
maps of the current as a function of Esd and ep. In one of 
the panels, we also report the difference between the cur- 
rent obtained within the quantum approach and that of the 
semi-classical calculation. This difference is higher where 
the Zener tunneling is more important. For metallic tubes, 
the current does not depend on the Fermi level because, 
within the linear band approximation, the electronic density 
of states is a constant. Note that for sufficiently high ep val- 
ues the semi-classical and the quantum results for the Qm 
tubes are almost indistinguishable and approach the results 
for the metallic tubes. 



C. Steady-state distribution functions 

It is instructive to compare the steady-state distribution 
functions of the electronic populations associated with the 
semi-classical or quantum approaches. Fig. [5] reports the 
distribution functions obtained for le = 50 nm and the 
Fermi level ep = 0.05 eV, that is for conditions in which 
the quantum effects are relevant. In Qm tubes, when the 
quantum approach is used both electron and hole contribute 
to the current, whereas in the semiclassical approach only 
electrons carry a current since the valence band is filled. In- 
terestingly, at high field the population of Qm tubes tends 
to that of metallic tubes only when the quantum approach 
is used. Indeed, in the metallic tube (which corresponds to 
the limit A — > 0) the tunneling between valence and con- 
duction band is total. The metallic limit can be reached 
only when one includes interband transitions. Fig. [5] also 
shows the imaginary part of the off-diagonal element of the 
density matrix Im(pi_i). This quantity is different from 
zero when quantum effects are relevant. One can remark 
that Zener tunneling turns on near k = Q (with a maximum 
at = ) and decays after few oscillations for fc > . 
The typical wave vector associated with these oscillations is 
cx 1/A, A being the electronic gap. 

VI. CONCLUSION 

We have presented a theoretical model based on the mas- 
ter equation to explore the quantum effects (inter-band 
transitions induced by electric field or Zener tunneling) on 
the transport properties of homogeneous quasi-metallic car- 
bon nanotubes (Qm). The presence and relevance of the 
Zener tunneling has been highlighted by comparing with a 
semi-classical Boltzmann transport approach, in which the 
inter-band transitions are neglected. We studied Qm tube 
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FIG. 6: (Color online) Current as a function of the electric 
field EsD and of the Fermi level ep for an elastic scattering 
le = 300 nm. For the quasi- metallic tube (Qm), calculations 
are done using the semi-classical or the quantum approaches. 
The difference between the semi-classical and quantum calcu- 
lations [dif. Sem-Quant (Qm)] is also shown and the current 
for this panel is multiplied by a factor 2. Calculations for the 
metallic tube are shown as a comparison. 
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FIG. 7: (Color online) Current as a function of Esd and ep 
for an elastic scattering le — 50 nm. Same caption as in Fig.[71 



with an electronic gap A = 50 meV in the presence of 
impurities or defects which influence the transport by pro- 
viding an elastic scattering channel for the carriers. Zener 
tunneling is relevant for small doping, when the Fermi en- 



FIG. 8: (Color online) Equilibrium distribution functions for 
conduction (pii) and valence bands (p_i_i) and imaginary 
part of the off-diagonal (Im(pi_i)) of the density matrix. Cal- 
culations are done with le — 50nm and ep = 0.05 eV, for two 
different electric fields Esd- The three panels in the left col- 
umn correspond to Esd = 0.5 V/^m (low filed), while those 
on the right column to Esd = 3.5 V/^m (high field). The 
two panels of the topmost row are calculations done for Qm 
tubes within the semi-classical (BTE) approach. The middle 
row panels are Qm tubes done within the quantum approach 
(ME). The bottom row panels correspond to metallic tubes. 
In this case p„„ (Ppp) is the distribution function of the band 
with negative (positive) velocity (see Fig. [!}. 



ergy lies in or close to the forbidden gap A. In absence of 
elastic scattering (in high quality samples), the small size 
of the band gap of Qm tubes enables Zener tunnelling for 
realistic values of the the electric field (above ^ 1 V//im). 
However, for such electric fields the scattering with the op- 
tical phonons (whose effect is underestimated in the present 
approach) tends to mask the Zener current. On the other 
hand, the presence of a strong elastic scattering (in low 
quality samples) further decreases the values of the field 
required to observe Zener tunnelling. Indeed, for elastic- 
scattering lengths of the order of 50 nm, Zener tunnelling 
affects the current/voltage characteristic already in the lin- 
ear regime. In other words, Zener tunneling is made visible 
by defects in Qm tubes. This result is similar to what has 
been recently shown for graphene^. 
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Calculations were done at IDRIS, Orsay, project 096128. 
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